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ABSTRACT 

We report on Chandra imaging observations of the Galactic Unidentified 7-ray source GEV J 1809-2327, com- 
paring the X-ray images with new VLA 1.46 GHz and 4.86 GHz maps. The X-ray images reveal a point source 
connected to a non-thermal X-ray/radio nebula, supporting a pulsar/wind model for the 7-ray emitter. We also 
detect numerous X-ray sources from the young stellar association in the adjacent HII region S32. 

Subject headings: gamma rays: observations - pulsars: general - supernova remnants 
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1. INTRODUCTION 

The EGRET instrument on the Compton Gamma Ray Obser- 
vatory detected many bright sources along the Galactic plane, 
most of which remain unidentified. A number of these sources 
have high significance in E >GeV photons (Lamb & Macomb 
1997), providing relatively good source localizations. The key 
to making progress on identifications is to search for lower en- 
ergy, particularly hard X-ray and radio, counterparts. Roberts, 
Romani, & Kawai (2001, RRK) have carried out an ASCA sur- 
vey of the bright GeV sources to search for counterparts and 
test the nature of the GeV emitters. 
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FIG. 1. — VLA 21cm map with ACIS 0.5-8keV smoothed contours. The 
ACIS-I FOV (square) and the brightest sources of Table 1 are also indicated. 

GeV J1809-2327(=2EG J181 1-2339) is one of the brightest, 
best-localized unidentified plane sources. It is one of a subset 
of sources showing evidence of 7-ray variability, with an ex- 
tended hard-spectrum X-ray counterpart argued to be a pulsar 



synchrotron wind nebula, or PWN (RRK 2001). This source is 
also interesting because it is adjacent to the dark cloud Lynds 
227 and to the HII region S32 with an embedded association of 
young high-mass stars. These are plausibly associated with the 
X-ray source and at ~ 1 .8kpc provide a possible birth site for a 
pulsar. Oka et al. (1999), in a study of the molecular gas in this 
region, found morphological support for this association and 
hypothesized that the 7-ray photons were produced when TeV 
pulsar electrons penetrated the dense molecular gas, suffering 
relativistic bremsstrahlung losses. 

Studies of this complex region were hampered by the modest 
angular resolution of the ASCA GIS (~ 3' half power diameter), 
so we have obtained high resolution ACIS images to separate 
out the point source contribution to the X-ray complex. This is 
abetted by new VLA continuum maps. These data support the 
pulsar/PWN hypothesis and give new constraints on the origin 
of the high energy emission. 

2. OBSERVATIONS AND DATA ANALYSIS 
2.1. CXO Imaging and Source Detection 

We observed the gamma-ray source GeV J 1809-2327 on 25 
August 2000 with the ACIS detector (Burke et al. 1997) on 
the Chandra X-Ray Observatory (Weisskopf et al. 2000) for 
~ 9.7ks. The exposure had the imager chips 10-3 at primary 
focus with additional off-axis coverage from chips S2 and S3. 
Integrations were made in TE mode. Basic analysis of the (re- 
processed) data was accomplished using CIAO v2.2.0.1. We 
found no evidence of background event rate flares; therefore, 
the entire dataset was used in the analysis. 

Our primary X-ray source list was established using wavde- 
tect with a significance threshold of 10~ 7 . Forty four sources 
are detected with a quoted significance > 3.0, and in Table 1 
we list those sources with a significance of greater than 8.0. 
The most significant of these are marked in Figure 1. Com- 
parison of our optical and X-ray frames allows us to boresight 
positions to ~ 0.3". Several of these sources have matches in 
SIMBAD and the USNO A2 catalogues. Many of these lie in 
the S32 stellar association to the south of our GeV counterpart 
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(Figure 1). An overlay of the X-ray image with POSS II scans 
reveal roughly a dozen additional matches, either with optical 
m > 20 or lower (albeit still high) X-ray significance. 

We tested these sources for variability, comparing the indi- 
vidual source photon arrival times to the total image counts 
(mostly background) arrival times, using the K-S statistic (Press 
et al. 1992). Sources with a variability probability greater than 
~ 95% we regard as likely variables, while most sources are 
consistent with a constant flux. 

As in earlier ASCA images, it is evident that in the northern 
half of the field there is also diffuse X-ray emission. After 
exposure correction and application of adaptive smoothing fil- 
ters we find a wedge of diffuse emission in our field with 
point source 5 at its apex. Figure 1 shows contours from the 
smoothed image overlaid on a 1.4GHz continuum map. At 
very low surface brightness below the faintest contours shown, 
the diffuse X-ray emission seems to extend to the lower left, 
which might place point source 5 towards the center, rather 
than apex of the diffuse flux. However, the brightest emission 
to the NW correlates well with the radio continuum flux. 
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FIG. 2. — DSS image around source 5. Smoothed ACIS contours show the 
point source and 30" 'jet' . The ellipse bounds our nebular extraction region. 

Figure 2 overlays X-ray contours on the POSS2 R image. 
There is no optical point source to the plate limit within ~ 1 .5" 
of the X-ray position. Optical extinction is modest at source 5, 
but increases rapidly ~ 30" to the NE, toward the edge of Lynds 
227. There is significant emission in a narrow jet or trail from 
point source 5 to the diffuse nebula. In the raw photon counts 
it shows as a narrow linear feature with a brighter knot ~ 15" 
from the point source. 

2.2. X-ray Spectral Fitting 

For the brightest sources (> 80 counts), we have attempted 
individual spectral analyses, creating PI spectra, Ancillary Re- 
sponse Files, and Response Matrices. Because of the limited 
counts and the limited ACIS spectral resolution, several spec- 
tral models were, in general, formally acceptable. We chose to 
fit the bright stellar sources in S32 with an absorbed Raymond- 
Smith plasma model. Fit parameters and multi-dimensional er- 
ror estimates are listed in Table 1 . 



Our earlier ASCA observations had longer exposure (80ks) 
and larger effective area. We thus compare present spectral es- 
timates with the higher statistical significance (but spatially un- 
resolved) estimates in RRK. The ASCA extraction aperture was 
quite large (4' radius), so background (particles and Galactic 
X-rays) dominates the ACIS counts in this region. Accordingly, 
we also define a smaller extraction region (a 68" x 38" ellipse 
plus exclusion of the point source) covering the brightest part 
of the nebula (Fig. 2). In the larger aperture, we attempt to 
account for the particle-induced background by subtracting the 
ACIS 'blank-sky' images. For the smaller aperture, this blank 
sky subtraction made no significant difference in the source pa- 
rameters; for this and point source extraction regions, we skip 
this step. In all cases we carefully chose background regions 
on the ACIS-I chips that were free of obvious emission and 
spanned a similar range of CCD rows to the source aperture. 
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2.2.1. Nebula Analysis 



FIG. 3. — Count spectrum and models for point source 5. The excess counts 
above 3 keV are highly significant and require a hard model component. 

ASCA measurements of the diffuse emission determined a 
power-law index of T = 2.09+q [7 and an absorption of Nh = 
1 - 77 -oil x 10 22 cm" 2 (90% errors). Our fits to the ACIS counts 
in the small extraction region produces values consistent with 
these (Table 1 - quoted errors are multi-dimensional from the 
sherpa command projection). In the larger region we obtain 
power-law fit parameters several sigma away from the ASCA 
values. We believe that, even with the blank sky subtraction, the 
dominant particle contribution is not adequately removed in this 
large aperture. Unlike the small aperture, we also find that the 
ACIS fit parameters are significantly affected by the region cho- 
sen for the local background subtraction. If, however, we fix the 
spectral parameters at the ASCA-determined values and fit the 
counts in the large aperture, we obtain 3.6 x 10" 12 erg cm~ 2 s _1 
for the ACIS data, in good agreement with the flux measured in 
(RRK 2001). We therefore consider that the CXO data are con- 
sistent with the ASCA fits, but adopt the latter spectral indices 
as fiducial for the diffuse emission. 



2.2.2. Point Source Analysis 

For source 5, our counterpart point source, we fit a sim- 
ple absorbed blackbody. While statistically acceptable, this 
model cannot produce the counts above 3keV (Figure 3), which 
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given our small ~ 2" extraction aperture are highly signifi- 
cant (< 0.01 background counts and < 1 pile-up count are 
expected). Interestingly, several X-ray emitting young pul- 
sars have hard power law components (presumably magneto- 
spheric) superposed on soft thermal spectra (e.g. Saito et al. 
1997; Halpern & Wang 1997; Pavlov et al. 2001), with a typical 
index T = 1.5. We therefore assume this V and fit a powerlaw 
amplitude from the E > 3keV source counts. We then fix this 
power law and re-fit the soft BB component. Parameters for 
the combined model are listed in Table 1 . (with the combined 
fit smallest by ~ 2x), a result of the low count statistics in the 
high energy bins. 

2.3. Radio Continuum Imaging 

Mosaiced continuum observations of the GeV J 1809-2327 
field at 1 .46 GHz and 4.86 GHz were performed with the VLA. 
Details of the radio observations and full results will be pre- 
sented elsewhere (Roberts, Gaensler, and Romani in prepara- 
tion). The 1.46 GHz field (Fig. 1, 18" x 15" beam) shows dif- 
fuse emission associated with S32 and the X-ray PWN; both 
are also well covered in the 4.86 GHz mosaic. 

The S32 region has an inverted, partly thermal spectrum, 
with spectral index a w 0.1. The morphology is hollow cen- 
tered, following the contours of the infrared HII region. In con- 
trast, the wedge of continuum to the north has a non-thermal 
spectrum, with an approximate spectral index of a R « -0.3 ± 
0.1 and significant polarization. The emission follows our dif- 
fuse X-ray trail extending ~ 6' to the NW, where the X-ray 
emission fades below detectability. 

3. INTERPRETATION AND CONCLUSIONS 

Comparing our new observations with the pulsar+PWN hy- 
pothesis, we note that source 5 composite spectrum is charac- 
teristic of a young pulsar, with a hard (presumably magneto- 



spheric) power law and underlying thermal emission. Given 
the high kT = 0.18keV, our fit thermal flux implies a small ef- 
fective area of A eff w 4 x 10 11 (<i/l-8kpc) 2 cm 2 . At < 3% of a 
neutron star's surface, this could represent a polar cap heated 
by precipitation of magnetospheric particles. Alternatively, at- 
mosphere effects (Romani 1987; Pavlov et al. 2001) can pro- 
duce a hard thermal tail, although such a large fit kT is diffi- 
cult to produce for reasonable pulsar ages. At a more plausi- 
ble kT e ff < 0.05keV, full surface emission from cooling would 
be unobservable; our data only limit such a component to 
<1.2xl0 6 K. 

Both the power law emission and the hot thermal com- 
ponent can thus measure magnetospheric activity, and have 
been phenomenologically related to the spin-down luminos- 
ity. Becker & Triimper (1997) find that for pulsars observ- 
able by ROSAT, the flux scaled as L x (0.l -2.4keV) w 10" 3 £; 
after correcting for absorption our inferred 0.1 -2.4keV flux 
of w 2 x 10~ 12 erg cm" 2 s _1 corresponds to a spin-down lumi- 
nosity E = 8 x 10 35 if 2 8 erg s" 1 at the fiducial source distance. 
Saito et al. (1997) find that the (pulsed) ASCA luminosity scales 
as L(2- lOkeV) w 10 34 £ 3 3 f erg s _1 . Using our 2-10keV un- 
absorbed point source flux, we infer from this relationship 
E = 1.5 x 10 36 c/[ 8 3 erg s" 1 . For dipole spindown we have E w 
10 38 (Z?i2T4)" 2 erg s" 1 ; these luminosities indicate typical spin 
parameters for a 7-ray emitting pulsar of B12T4 « 10, where the 
surface dipole field is 10 12 Bi2G and the pulsar age is 10 4 74y. 

The broadband spectral energy distribution (SED, Figure 4) 
of our PWN candidate provides some useful constraints on its 
nature, if we interpret the radio/X-ray spectrum as synchrotron 
emission from a cooling electron population: 6N e /dT e = AT~ S 
with As = 1 at a break T B . This corresponding to a photon break 
frequency v B = 4.2 x \Q n B w Y\, with a typical PWN field of 
lOBio^G and a electron break at a Lorentz factor T B = 10 6 T 6 . 
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Our SED allows a modest range for the power law index s = 
2.85 ± 0.15. The data require a very tight correlation between 
s and the break frequency, Log(V B ) = 9.1 s— 13.3; however this 
allows a large range of break energies from the microwave band 
through the near IR. 

The conventional PWN picture identifies this break with the 
cooling of the electrons at r 4 = ^B^T^ 1 , associating this with 
the pulsar characteristic age. In this approximation (a homoge- 
nous, uniform field PWN) our data then can be used to elimi- 
nate Tb in favor of s, giving a PWN field estimate of 
Log(B 10 ) « 0.73-2/3Log(T 4 )-3.03(>-3). 

The total energy of the electron population, w lO^B^^erg, 
gives an upper limit on the initial pulsar period of P t < 

0.4Bjq +1)//4 s. Of course, adiabatic losses in the PWN flow likely 
decrease the present electron energy and require smaller P t . 
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FIG. 4. — Our PWN candidate's broad-band SED, with X-ray error regions 
for the full nebula (ASCA) and bright CXO core. The curve is the best-fit 
i = 2.85 cooling synchrotron model. 

We can check consistency with another crude estimate of 
the nebula field. Oka et al. (1999) argued that the 'fit' of 
the PWN candidate with Lynds 227 suggests confinement and 
pressure equilibrium. The CO line-width inferred pressure of 
~ 1.4 x 10" n g/cm/s 2 implies a nebular equipartition field of 
~ 20/iG. However, since Lynds 227 brackets only on one side, 
the wind is likely only partly confined and B may be larger. For 
our best fit s, the equipartition B suggests r 4 <~ 20, a Geminga- 
like pulsar. Our present poor constraint on v B allows ages 
~ 3 x 10 4 - 10 6 y, including B\2T\ « 10 as estimated above. If 
the PWN B exceeds 20fiG, even smaller r are permitted. 



We can also follow the pulsar surface field Byi out to the ter- 
mination shock in a simple 1-D picture. The magnetic field 
at r s = 9 s d, the wind shock radius, is B s ~ 3B*rl/(rj jC r s ) ~ 
130(Bi2T4<fi.86'')~ 1 AiC" Beyond r s the field evolution depends 
on the poorly understood wind magnetization, but assuming 
B s « 20^G constant in the postshock wind, for B^ta ~ 10 
we expect termination at 6 S s=s 40". Our ACIS image does not 
show a spherical wind with a subluminous zone at this 9 S , but 
instead suggests a jet or pulsar trail, with the ~ 15" knot in 
the jet/trail plausibly a termination shock. The overall dif- 
fuse X-ray morphology suggests a pulsar trail leading back to 
a PWN, as seen for PSR B 1757-24 (Kaspi et al. 2001), im- 
plying a birthsite to the NW. However, the nebula might also 
be powered through a jet, as for the asymmetric PSR 1509-58 
PWN (c.f. Gaensler et al. 2001). In this case the faint extended 
emission to the SE (Figure 1) may be the counterlobe. Un- 
der this interpretation, the S32 cluster at r « 5<ii.8pc provides a 
plausible pulsar birthsite, requiring a pulsar transverse velocity 
vj_ w 500<ii.8/T4km/s. Much higher S/N imaging or a proper 
motion measurement are needed to distinguish these possibili- 
ties. 

The 7-ray flux is clearly an additional spectral component 
(Figure 4), but its origin is puzzling. 3EG J 1809-2328 (= GeV 
J 1809-2327) is apparently one of the more variable Galactic 
plane sources (Tompkins 1999). Oka et al. (1999) suggested 
that the PWN-generated electrons penetrate the molecular 
cloud, generating GeV photons via relativistic Bremsstrahlung. 
As the PWN is well separated from the molecular gas at ~ 
(3')d ~ 1.6pc, one might expect variability no faster than the 
3//c~ 15y flow-crossing time. This is substantially longer than 
the EGRET variability timescale, suggesting an origin closer to 
the pulsar. Inverse Compton emission from the PWN termi- 
nation shock seems attractive, but both Galactic IR and local 
synchrotron emission fluxes fail by several orders of magnitude 
to account for the required soft photon energy density. Even 
the optical/UV emission from the OB stars in S32 contribute 
less than 10 5 5 L , failing to produce the required target photon 
density at the ~ 5pc PWN distance by > 300 x . The 7-ray flux 
is quite plausible for magnetospheric (curvature) pulsar emis- 
sion at the inferred E and distance, but variability is unexpected 
from an isolated pulsar. Confirmation or exclusion of the 7-ray 
variability, together with detection or limits on 7-ray pulsations 
is the best way to address this issue. Unless an X-ray period 
can be measured, this will not be resolved until until AGILE or 
GLAST make new sensitive 7-ray observations. 

This work was supported in part by CXO grant GOO- 1125 
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referee for a detailed critique. 



REFERENCES 



Becker, W. & Trumper J. 1997, A&A 326, 682 

Burke, B.E., Gregory, J., Bautz, M.W, Prigozhin, G.Y., Kissel, S.E., Kosicki, 
B.N., Loomis, A.H., & Young, DJ. 1997, Transac. Elec. Devices, 44, 1633 
Halpern, J. P. & Wang, F. Y.-H. 1997, ApJ, 477, 905 

Gaensler, B.M., Arons, J., Kaspi, V.M., Pivavaroff, M.J., Kawai, N., & Tamura, 

K. 2001, ApJ, submitted 
Kaspi, V.M., Gotthelf, E.V., Gaensler, B.M. & Lyutikov, M. 2001, ApJ, in 

press. 

Lamb, R. C. & Macomb, D. J. 1997, ApJ, 488, 872 

Oka, T., Kawai, N, Naito, T., Horiuchi, T., Namiki, M., Saito, Y., Romani, R. 

W., & Kifune, T. 1999, ApJ, 526, 764 
Pavlov, G.G., Zavlin, V.E., Sanwal, D., Burwitz, V. & Garmire, G.P 2001, 

ApJ, 552, L129 



Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, 
Numerical Recipes in Fortran, (Cambridge University Press:Cambridge) 
Roberts, M. S. E., Gaensler, & Romani, R. W. 2001, ApJ, (in press) 
Roberts, M.S.E., Romani, R.W., & Kawai, N. 2001, ApJS, 133, 451 
Romani, R.W. 1987, ApJ, 313, 718. 

Saito, Y, Kawai, N, Kamae, T., & Shibata, S. 1997, Neutron Stars and Pulsars: 
Thirty Years after Discovery, ed. N. Shibazzaki et al. , Tokyo, Universal 
Academy Press, 295 

Tompkins, W. 1999, PhD Thesis, Stanford University 

Weisskopf, M.C., Tananbaum, H.D., Van Speybroeck, L.P, & O'Dell, S.L. 
2000, Proc. SPIE, 4012, 2 



